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In our recent publication (D. O'Dell, et al, Phys. Rev. Lett. 
84, 5687 (2000)) we proposed a scheme for electromagneti- 
cally generating a self-bound Bose-Einstein condensate with 
1/r attractive interactions: the analog of a Bose star. Here 
we focus upon the conditions neccessary to observe the tran- 
sition from external trapping to self-binding. This transition 
becomes manifest in a sharp reduction of the condensate ra- 
dius and its dependence on the laser intensity rather that the 
trap potential. 



I. INTRODUCTION 

We have recently proposed 0] a scheme for induc- 
ing a 1/r gravitational- like attractive interatomic poten- 
tial in an atomic Bose-Einstein condensate (BEC) Q 
contained in the near-zone volume of intersecting tri- 
ads of orthogonal laser beams. For sufficiently strong 
self- "gravitation" the BEC becomes self-bound. In this 
unique, novel regime the 1/r attraction balances the out- 
ward pressure due to the zero point kinetic energy and 
the short range s-wave scattering. Here we focus upon the 
transition from external trapping to self-binding. This 
transition becomes manifest in a sharp reduction of the 
condensate radius and its dependence on the laser in- 
tensity rather than the trap potential. We analyze the 
conditions for the observability of the self-binding tran- 
sition: the threshold laser intensity (Sec. |n|), the bounds 
on the num ber of atoms imposed by the near-zone con- 
dition (Sec. Ill), as well as the loss rates (Sec. IV). Sec. 
^ summarizes the findings. 

II. SELF-BINDING THRESHOLD INTENSITY 

A. Threshold condition 

We need to find a situation where the mean-field self- 
"gravitation" energy associated with the near-zone laser- 
induced attractive 1/r potential can become (at least) 
comparable with the short-range s-wave scattering en- 
ergy. To this end, we examine the mean-field solution for 
a condensate of atoms interacting via Thirunamachan- 
dran's isotropic two-atom potential J3|, obtained by di- 
rectional averaging of the laser-induced dipole-dipole po- 
tential. This potential has the form 
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where f 
Al, and 



t/Al, is normalized to the laser wavelength 
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/ being the sum of the intensities of all the lasers, and a 
the atomic polarizability. The potential begins to oscil- 
late (i.e. becomes alternatingly repulsive and attractive) 
at distances beyond ~ 0.36Al. However, this potential 
can support a self-bound condensate with a larger radius 
as shown below. 

We use the mean- field approximation (MFA), as em- 
bodied in the following generalized Gross-Pitaevskii 
equation Q , to calculate the ground-state order parame- 
ter 'J'(R) of a BEC subject to a laser-induced interatomic 
interaction 
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where m is the atomic mass, V cxt (R) = mui^R 2 /2 is 
an isotropic external trap potential (which will be con- 
sidered negligible — see below), and 14 C (R) is the self- 
consistent potential 



V SC (R) = gp(K) + / d 3 R' L7 iso (R - R) p(R') 
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where p(R) = ^ 2 (R) is the density and g = Airaft 2 /m, a 
being the s-wave scattering length. 

In cold dilute atomic BECs with short-range s-wave 
scattering, the validity of the MFA (i.e. the Gross- 
Pitaevskii equation [||), is well established providing 
pa 3 <C 1. However, the MFA is also valid for the long 
range repulsive coulomb potential, +u/r, provided many 
atoms lie within an interaction sphere with a Bohr-type 
radius, a* = h 2 /mil, so that pa 3 > 1 This condition 
means that the potential must be weak. Remarkably, 
sclf-gravitating BECs simultaneously satisfy both of these 
MFA validity conditions as can be readily verified using 
the ensuing expressions. 

The condensate radius can be studied using the vari- 
ational wavefunction ^> W (R) = yiV ey^){— R 2 / 2w 2 X 2 /) / 
(ttw 2 X 2 /)* , where w is a dimensionless variational param- 
eter giving the width of the condensate. The variational 
solution in the limit of negligible kinetic energy (Thomas- 
Fermi limit) yields a self-bound condensate, i.e. finite w 
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(see Fig. 1 and Fig. 2 below), if the laser intensity exceeds 
the following threshold value (in S.I. units) 
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Here Iq is the total intensity supplied by all the laser 
beams: for a triad each laser should have 1/3 of the 
above value and for the 6 triad configuration (l) 12 of 
the lasers should have 1/15, and the remaining 6 should 
have 1/30, of the above value. The threshold Iq signifies 
the equality of the gravitational-like potential and the 
s-wave scattering potential. 

With an intensity 1.5 times the threshold value (Eq. 
(0)) (arrow in Fig 2) the expectation value of the rms 
condensate radius R rms = sj (R 2 ) is a fraction of the 
laser wavelength Al {Rrms ~ 0.43 x Al). The condensate 
is less and less confined as one approaches the thresh- 
old (||) — see Fig. 2, from above. Increasing the intensity 
/ reduces the condensate radius which becomes, in the 
asymptotic limit, proportional to 1/vT. Thus the de- 
pendence R rms ~ (Jo/J)^ 2 Al is a distinct experimental 
signature of self-binding. 




FIG. 1. (a) Variational mean field energies per particles 
in the case of negligible kinetic energy (TF-G regime) and 
Xh/Na <C 1 plotted versus the trial size w for different values 
of I /Iq. (b) Equilibrium value of w versus I /Io in the limit 
of negligible kinetic energy (Thomas- Fermi limit). Only for 
J/Jo > 1 are self-bound variational solutions (having mini- 
mum at finite w) observed. Inset - Schematic phase portrait 
of the transition from unbound to self-bound regime for neg- 
ligible external trapping is plotted versus log 10 (Xi J /Na) and 
log 10 (/// ). 

At the threshold intensity an external harmonic trap 
becomes negligible when plgX^a 1, where I 2 = h/muio 
and p is the density. As the laser intensity is in- 
creased beyond this value the trap becomes increasingly 
"irrelevant" — it is not necessary to turn it off to access 
the TF-G regime, where r _1 and s-wave scattering dom- 
inate. 

The threshold Iq (||) is evaluated neglecting the ki- 
netic energy. The role of kinetic energy can be discussed 
in terms of X^/Na (approximately the ratio between 
the kinetic energy Nh 2 /mX^ and the scattering energy 



N 2 h a/mA L ), as shown schematically in the phase por- 
trait in Fig. 1 (drawn for negligible external trapping) 
which can modify the threshold for self-binding. The G 
regime, representing the purely "gravitational" counter- 
part of the TF-G regime, where only "self-gravitation" 
and kinetic energy play a role jij (as in a Bose star ||) 
is accessed when 
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that implies 1 < X^/Na. 

At this point the variety of choices can be mainly di- 
vided into two categories: i) to work with long laser wave- 
lengths in order to contain many atoms within the near 
zone, at the price of very high threshold power; ii) use 
laser wavelengths moderately detuned from an atomic 
resonance, so as to benefit from the increased polarizabil- 
ity, at the price of considerably fewer self-bound atoms. 



B. Long- wavelength (static polarizability) threshold 

The threshold intensity (Eq. (^J)) is independent of the 
laser wavelength Al, as long as the dynamic polarizabil- 
ity a(q) is too. The io threshold takes the following zero- 
frequency (static) values: Iq = 5.65 x 10 9 Watts/cm for 
sodium, Iq = 8.19 x 10 8 Watts/cm for rubidium. It is 
sufficient to use 20 W x 3 beams of Nd:Yag lasers focused 
down to 10 /im for rubidium to exceed the threshold. By 
contrast, we require multi kW CO2 lasers focussed down 
to 100 /im for the same purpose. A laser beam with a 
gaussian profile focussed to 10 Al would exert a large in- 
ward radial dipole force on each atom, so non-gaussian 
optics giving a very flat intensity profile || over the con- 
densate region may be required in the long-wavelength 
(static) case. There remains the problem of random noise 
in the intensity profile, but fortunately this can only ex- 
ist on scales larger than the wavelength and so may be 
overcome. 

An additional option is to reduce the scattering length 
a (to which the threshold intensity (|^) is proportional). 
This is possible in the vicinity of (but somewhat off) a 
Feshbach resonance, as demonstrated experimentally [0: 
reduction of a, and correspondly Iq, by one to two orders 
of magnitude would eliminate the need for non-gaussian 
optics in the static polarizability case. 



C. Moderate-detuning threshold 

Using a moderate detuning from an atomic resonance 
one can increase the polarizability by many orders of 
magnitude compared to its zero frequency value. 

In a recent experiment on superradiance || the laser 
was red detuned by 1.7 GHz from the 3S]/2, F=l — > 
3P3/2, F = 0,1,2, transition of sodium. With this de- 
tuning, the polarizability in cgs units is a — 3.534 x 
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10~ 18 cm 3 , which is ~ 1.5 x 10 5 times the static value 
of the polarizability. The threshold intensity (^) is then 
reduced by a factor w 2.3 x 10 10 compared to the static 
polarizability case, becoming 1$ 262 mW/cm 2 for 
sodium, which is close to the values used in Ref. M . 

With this value of threshold intensity the gradient 
forces can be negligible if the focal spots of the lasers 
are much wider than Al- 



D. Moderate-detuning saturation and repulsion 

The potential (Q) is the result of a 4th order, two- 
atom, QED process |J, valid when the laser is far de- 
tuned from any atomic transitions. This means that the 
initial absorption of a laser photon and the subsequent 
intermediate steps are virtual processes (which are most 
significant in the near-zone), followed by photon emission 
back into the original laser mode. A different process can 
take place when the laser is on-resonance. Genuine ab- 
sorption of a laser photon by a single atom (measured 
by the saturation), followed by spontaneous emission of 
this real photon is a process that radiates energy. If an- 
other atom absorbs this radiation then in the far-zone it 
feels a repulsive Coulomb-like force F repu i s = K/r 2 [ pr| , 
which has been recently measured in rubidium molasses 
p2[ . For moderate detuning, can this force counteract 
our attractive gravitation- like force -Fgrav = —u/r 2 7 

For detuning 5 much larger than both the Rabi fre- 
quency Q and the linewhidth 7 of the resonance, the 
saturation parameter s — Id 2 /(each 2 S 2 ) jlQ], where d 
is the dipole matrix element, becomes independent of the 
detuning when calculated at the threshold intensity @ 
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It is then found that ]l2j K w cr^/ s ft 4 /(16c<5 2 ), where 
is the resonant absorption cross section and I s is the 
corresponding saturation intensity. On comparing this 
expression with u (Eq. @), we find that, in terms of the 
saturation parameter s, 
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For the sodium transition and 1.7 GHz detuning referred 
above, Eq. ([?]) yields very small value s » 0.0003. This 
implies that under the moderate-detuning conditions dis- 
cussed above, the repulsive force has a negligible effect on 
self-binding. 



have the maximum density of some 10 15 atoms/cm 3 . Us- 
ing the gaussian wavefunction one can have of the order 
of 10 6 or 10 3 atoms in the condensate irradiated by a 
CO2 laser or Nd:Yag laser, respectively (see Fig. 2). 

The price of moderately detuned wavelengths (~ 
.589/xm for sodium) is the small number of atoms in- 
volved. With an intensity I » 1.5/o the atom cloud con- 
tains w 40 atoms as the peak density ranges from 10 15 
to 10 16 atoms/cm 3 . Although this number is small, it is 
sufficient to demonstrate the self-binding effect. 

For given values of /, a, a, and m, we are either in the 
G regime or the TF-G regime, depending on whether the 
number of atoms N is smaller or larger than the number 

Abordcr ~ "J 3nh~ 2 / (2mua) which corresponds to the 
line separating the two regions in the inset of Fig. 1. 

It so happens that 40, the lower estimate of the number 
of self-bound sodium atoms obtainable in the moderate- 
detuning regime, is very close to A^order- This is an in- 
teresting region, because both the kinetic energy and the 
s-wave scattering are significant and together with the 
r _1 attraction determine the condensate properties. 
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FIG. 2. Range of numbers TV of Na condensate atoms as a 
function of Al that are compatible with a TF-G or G solution. 
The density is 10 15 -10 16 atoms/cm 3 and the intensity is 1.5 
times the threshold intensity (^). The region above 10 16 cm~ 3 
corresponds to excessive density. The vertical long-dashed 
line corresponds to the moderate-detuning choice discussed 
for Na. 



IV. LOSS RATES 



III. NUMBER OF SELF-BOUND ATOMS 

A key experimental restriction on self-binding is that 
the atoms should be in the near-zone to feel the 1/r po- 
tential: a condensate smaller than the laser wavelength 
limits the number of atoms involved. Let us assume we 



A. Spontaneous Rayleigh losses 

The single-atom Rayleigh scattering rate r ray leads 
to depletion of the condensate. The probability ampli- 
tude for inelastic scattering from the ground state |0) of 
the near-zone condensate to any excited state \n) due 
to an external field with wavevector q is proportional to 
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VN J2n^o( n \( c l ' r )|0)- Hence, for sample sizes less than 
a wavelength we expect the spontaneous Rayleigh scat- 
tering rate to be reduced by a factor at least as small 
as (qR rms ) 2 , analogously to the Lamb-Dicke effect [[l|. 
The lifetime of the condensate, when determined from 
spontaneous Rayleigh scattering alone, is estimated to 
be 



B. Intereference losses 
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Since r ray = Iq 3 a 2 /(3/iEqc) [||, it can be expressed in 
terms of the electromagnetically induced energy U (r) — 
—u/r of a single pair of atoms separated by a distance 
equal to the wavelength 
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where u is defined in Eq. (|2|). Using this relation, we can 
compare the upper bound on the condensate lifetime set 
by Rayleigh scattering with the time scale of the dynam- 
ics, the requirement being that the system exists long 
enough to equilibriate. In the TF-G and G (self-bound) 
regions a characteristic time scale for the dynamics is 
provided by the following "plasma" frequency 
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where /O p0 ak is the peak density. We can express lu p in 
terms of the recoil energy Er = h 2 q 2 /2m (q being the 
mean laser wavelength) and the Rayleigh scattering rate 
r ray using Eq. <M) 



where the factor 
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is asymptotically equal to 1 in the G region and 
N/N boldeI in the TF-G region. It follows from (|l2|) that 
the characteristic oscillation frequency lu p can be much 
bigger than r ray , by a factor proportional to N 2 or N 1 / 2 
in the G or TF-G region, respectively. Thus the lifetime 
can be considerably longer than the characteristic time 
scale of the dynamics. 

Even for the small number of 40 sodium atoms in the 
self-bound moderate-detuning regime (I = 1.5 X Iq, S = 
1.7 GHz), for which the recoil energy is Er/Ti — 1.57x 10 5 
s~ x and r ray = 1.58 x 10 4 s _1 , we find lu p w 20 x r ray . 
This implies that several oscillation periods of the self- 
bound condensate can occur within the Rayleigh lifetime. 



We revisit the expressions for the loss rate T^torf due 
to multi-beam interference as obtained in Q . We can ex- 
press Tinterf in terms of the recoil energy Er and Rayleigh 
scattering rate r ray as in Sec. IV A 
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where Q is the relative det uning of beams in the triad. 
In the example given in Sec. IV A above, T^tcrf turns out 
to be few times bigger than r ray , when f2 is chosen to be 
of the order of uj„. 



V. CONCLUSIONS 

Our main conclusion is that at least the TF-G self- 
bound region is experimentally accessible, although such 
an experiment would be challenging. Moderate detuning 
is preferable to the longer wavelength case due to the 
huge enhancement in the polarizability, but it allows the 
self-binding of few (less that 100) atoms. If the scattering 
length were reduced via a Feshbach resonance then this 
would further facilitate the self-trapping of many more 
atoms using near-infrared lasers. 
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